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O ■ ABSTRACT 

^ ' We present a high-resolution (~5") image of the nucleus of M82 showing the 

presence of widespread emission of the formyl radical (HCO). The HCO map, 
^ ' the first obtained in an external galaxy, reveals the existence of a structured 

I disk of ~650pc full diameter. The HCO distribution in the plane mimics the 

ring morphology displayed by other molecular /ionized gas tracers in M82. More 
precisely, rings traced by HCO, CO and HII regions are nested, with the HCO 
ring lying in the outer edge of the molecular torus. Observations of HCO in 
galactic clouds indicate that the abundance of HCO is strongly enhanced in 
the interfaces between the ionized and molecular gas. The surprisingly high 
overall abundance of HCO measured in M82 (X(HC0)~4 10~^°) indicates that 
its nuclear disk can be viewed as a giant Photon Dominated Region (FDR) of 
~650pc size. The existence of various nested gas rings, with the highest HCO 
abundance occurring at the outer ring (X(HCO)~0.8 10~^), suggests that PDR 
chemistry is propagating in the disk. We discuss the inferred large abundances of 
HCO in M82 in the context of a starburst evolutionary scenario, picturing the 
M82 nucleus as an evolved starburst. 
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1. Introduction 

M 82 is the closest galaxy experiencing a massive star formation episode (Rieke et al. 
1980; Wills ct al. 1999). Its nuclear starburst, located in the central Ikpc, has been the 
subject of continuum and line observations made in virtually all wavelengths from X-rays to 
the radio domain. These studies indicate that the high rate of supernova explosions and the 
strong UV radiation fields have heavily infiuenced the physical properties and kinematics of 
the interstellar medium in M82. The high supernova rate has created a biconical outfiow 
of hot gas (Bregman, Schulman, & Tomisaka 1995; Shopbell & Bland-Hawthorn 1998) also 
observed in the cold gas and dust (Alton, Davis, & Bianchi 1999; Seaquist & Clark 2001). 
The discovery of a ~500 pc molecular gas chimney and a giant supershell in M 82, detected in 
SiO, indicates the occurrence of large-scale shocks in the disk-halo interface of the starburst 
(Garcia-Burillo et al. 2001). Furthermore, there are evidences that the strong UV-fields have 
created a particular physical environment in the molecular gas reservoir of M 82 (Stutzki et 
al. 1997; Mao et al. 2000; Weiss et al. 2001). The common picture emerging from these studies 
is that the bulk of CO emission in the nuclear disk of M 82 comes from moderately dense 
(with n(H2)'^10^-10^cm~^) Photon Dominated Regions (PDR). However, these conclusions 
are model-dependent and not free from internal inconsistencies (Mao et al. 2000). 

Observational evidence supports that the emission of the formyl radical (HCO) mainly 
arises from PDR at the interfaces between the ionized and the molecular gas in our Galaxy. 
After the first detection of HCO by Snyder, HoUis, & Ulich (1976) ulterior searches have 
confirmed that this radical is associated with regions where chemistry is driven by an en- 
hanced UV radiation field (HoUis & Churchwell 1983; Snyder, Schenewerk, & HoUis 1985; 
Schenewerk et al.l988; Schilke et al. 2001). On the theoretical side, several chemical models 
have been put forward to account for the observed abundances of HCO in PDR (de Jong, 
Boland, & Dalgarno 1980; Leung, Herbst, & Huebner 1984; Schilke et al. 2001). 

We have chosen HCO, a privileged tracer of PDR chemistry, to investigate the infiuence 

of the M82 starburst on its molecular gas reservoir. Sage & Ziurys (1995) tentatively detected 
the emission of HCO at 3mm in M 82, using single-dish observations. However, the HCO lines 



^Bascd on observations carried out with the IRAM Plateau de Bure Interferometer. IRAM is supported 
by INSU/CNRS (France), MPG (Germany) and IGN (Spain). 
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appear blended with SiO and H^^CO+ lines in their spectrum. Furthermore, the low spatial 
resolution (50" ) of their single-pointed map did not allow to infer the spatial distribution of 
HCO. In this Letter we present a high-resolution (~5") image of the emission of HCO in 
the nucleus of M 82. The interferometer HCO map of M 82, the first obtained in an external 
galaxy, shows unambiguous evidence that the whole nuclear disk has become a giant PDR of 
~650pc size with a total HCO abundance of ~410~^°. The enhancement of HCO presents 
spatial variations which depend on the distance to the most prominent HH regions of the 
M82's starburst. 



2. Observations 

The HCO observations of M 82 were carried out with the IRAM array at Plateau 
de Bure (France) during 1999 June, using the 4-antennas CD set of configurations. We 
observed four of the strongest 'NK-K+=^orOoo hyperfine lines of HCO. These transitions 
form a 3mm quartet: J=3/2-l/2, F=2-l(86.671 GHz) and F=l-0(86.708GHz); J=l/2- 
1/2, F=l-1 (86.777 GHz) and F=0-l(86.806GHz). We observed simultaneously the J=2-l 
line of SiO(86.847GHz), and the J=l-0 line of Hi3C0+ (86.754 GHz). The one-field ob- 
servations are made centering the 55" primary beam of the array at q;j2ooo=09'^55™51.9* 
and 5j2ooo=69°40'47.1", which corresponds to the 2.2/1 peak (Joy, Lester, & Harvey 1987). 
Calibration of visibilities was performed as explained in Garcia-Burillo et al. (2001). The 
synthesized beam is almost circular (5.9" x 5.6", PA=105°). The rms noise level in 2.5 MHz 
wide channel maps, derived after subtraction of the continuum emission, is ImJy/beam 
(5mK). We assume a distance to M82 of D=3.9Mpc (Sakai & Madore 1999); the latter 
implies 1" ~20 pc. We define the x and y axes to run parallel to the galaxy major and minor 
axes, respectively, where x >0 eastward and y >0 northward. The major axis has PA=70°; 
(x,y) offsets are referred to the 2.2/i peak. All through the paper, line intensities are not 
corrected from primary beam attenuation. Note however that the bulk of the discussion is 
based on the study of the ratios of lines which are equally affected by this bias. If they were 
to be corrected, HCO and H^^CO"'' intensities should be both multiplied by a factor ~1.3 at 
the edge of the M82's nuclear disk (at a~17" radius). 



3. The HCO map 

Figure la represents the velocity-integrated intensity map obtained for the strongest 
hyperfine component (F=2-l) of the J=3/2-l/2 line of HCO. This image shows that the 
emission of HCO is widespread in the nuclear starburst of M82. The bulk of the emission 
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comes from a highly-structured disk of ~650pc full diameter and barely resolved vertical 
thickness. The HCO distribution in the plane mimics the ring morphology displayed by 
other molecular gas tracers in M82 (Mao et al. 2000). Two emission peaks at (x,y) = (+16" ,- 
2") and (-15", 0") locate the eastern and western lobes of the HCO ring viewed edge-on. 
A central clump at (+2",0") is close to the 2.2fi peak identified as the dynamical center 
(Garcia-Burillo et al. 2001). With lower statistical significance, we detect out-of-the plane 
emission in a clump at (-9", 8"). However, the full-size of this structure is smaller than the 
beam, and therefore its reality is doubtful. 

Figure lb shows the spectra of HCO, SiO and H^'^CO^ integrated in the regions of 
the two ring maxima (see Figure la). Gaussian profiles arc fitted for the identified lines. 
Together with the strongest F=2-l fine, emission of the F=l-0 component is detected at 4(7 
level in the two ring lobes. We have also detected the F=l-1 and F=0-1 transitions in the 
western and eastern lobes, respectively. We estimate an average F=2-l/F=l-0 fine ratio of 
~3±1 (~1.5±0.5) for the western (eastern) lobe. Within the errors, these ratios are close to 
the theoretical line ratio of ~1.7, expected for optically thin emission. 

Figures 2a-b show the HC0(F=2-1) integrated intensity contours overlaid with the 
zero moment maps of H^3C0+ (1-0) (Garcia-Burillo et al. 2001) and ^^00(2-1) (Weiss 
et al. 2001). Although all images reveal a ring-like distribution of molecular gas in the 
nucleus of M82, the separation between the E/W lobes differ significantly between HCO 
(~32") and both ^^CO and Hi''^CO+ (~25"). The HCO emission extends farther out in the 
disk, especially in the eastern lobe. The HCO torus encircles the molecular torus traced by 
H^^CO^ and CO. Most remarkably, the ring of HII regions (of 10"-15" size) identified by 
Br7 (Larkin et al. 1994), Nell (Achtermann & Lacy 1995), and radio- recombination lines 
(e.g. the H41q; map of Seaquist et al. 1996) hes just inside the CO molecular ring. The 
HCO, CO and HII rings are nested with HCO lying in the outer edge of the disk. The 
central HCO maximum is also displaced from the ^^CO/H^^CO+ hot spot at (-3",0"). At 
smaller scales also, the peaks of HCO appear to avoid the brightest HII regions (see section 
4). To interpret the described morphology we will first estimate the spatial distribution of 
the abundance of HCO in M 82 and discuss different scenarios accounting for it. 

Kinematics of HCO gas in M82's disk show no significant departures from the general 
rotation pattern typically displayed by other molecular gas tracers in this galaxy. 
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4. The HCO abundance map in M 82 



Figure 3 represents in gray scale the HCO(F=2-1)/H^'^CO+(1-0) intensity ratio map 
in the nucleus of M82, derived using a Qa (Aa) clipping for Hi3CO+(l-0) (HC0(F=2-1)). 
The inferred ratio, accurate to 30% varies between ~0.50±0.15 (at the outer edges of the 
ring, and in the central peak) to ~0.15±0.04 (in the region between the HCO maxima), with 
an average ratio of ~0.25±0.07 for the whole disk. In order to derive the HCO to H^^CO"^ 
column density ratio we need to make several assumptions about both the optical thickness 
and excitation temperature of the lines. Based on studies of HCO emission in galactic 
PDR (Schenewerk et al. 1986, 1988) it is plausible to suppose that the HCO lines should 
be optically thin also in M82. The hyperfine line ratios measured in M82 seem to confirm 
these expectations. As excitation temperature for HCO we will assume Tcx=10-ft'; this value 
was derived from a multitransition analysis made by Snyder et al. (1985) on NGC2024. 

For H^'^CO^ we will also consider optically thin emission (Garcfa-Burillo et al. 2000) 
and the same excitation temperature as assumed for HCO. These are reasonable guesses, 
especially for Tex, as the two molecules have similar critical densities for the examined transi- 
tions. In this case the calculation of HCO to H^^CO+ column density ratio is straightforward 
using the expression (Schenewerk et al. 1988): 



where N is the total column density, gjj is the degeneracy of the upper hyperfine level 
{gu=5 for F=2), I is the integrated intensity, and A is the Einstein coefficient of the tran- 
sition. We derive values for N[HCO)/N{H^^CO^) that range from 2 to 8, with an average 
value of 3.6 for the whole disk. Adopting an average fractional abundance for H^^CO"*" of 
~10-i° (Garcia-Burillo et al. 2000, 2001), we derive X(HC0)~4 10"^° in M82. 

Schenewerk et al. (1988) measured a value of N{HCO)/N{H^^CO+) -9.7 in the HII 
region NGC2024. On the contrary, values oi N{HCO) / N{H^^CO^) significantly lower than 
1, i.e. an order magnitude below those found in NGC 2024, have been reported by Schenewerk 
et al. (1988) in galactic clouds without developed HII regions or having no indication of star 
formation. More recently, Schilke et al. (2001) have searched for HCO in a reduced sample 
of prototypical PDR. The estimated N{HCO)/N[H^^CO^) abundance ratios range from 
~30 (in the Orion bar) to ~3 (in NGC 7023). The largest HCO abundances are found in 
the Orion Bar: the paradigm of interaction between an HII region (M42) and its parent 
molecular cloud. Comparing our global estimates of N{HCO)/N{H^^CO~^) ~3.6 in M82 
with the values measured for the prototypical PDR in our Galaxy, we conclude that the 
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nuclear disk of M 82 can be viewed as a giant PDR. Further analysis of the spatial variation 
of the N(HCO)/N[H^^CO~^) ratio shows that the chemistry driven by the strong UV fields 
has produced the highest HCO enhancement in the outer edges of the starburst ring; the 
central clump also shows a similar HCO enrichment {N{HCO)/N{H^^CO^) ~7-8). These 
regions, with HCO abundances approaching those of NGC 2024, lie close to the loci of HII 
regions in M 82 but do not quite coincide spatially with them. A mere inspection of Figure 3 
indicates that the HCO-enriched molecular clouds surround the ionized gas as shown by the 
strongest [NeH] emission: X(HCO)/[Nen] peaks avoid each other and are shifted by offsets of 
~50-150pc. This might seem paradoxical as a close correlation between HE regions and PDR 
could be expected at these scales. As discussed below, this distribution could be explained 
however as the final result of propagation of PDR chemistry in M82. 



5. Discussion 

The enhancement of HCO in PDR is mostly an observational fact, reproduced with 
uneven success by models. PDR models published so far have tried to account for the 
derived abundance of HCO using either gas-phase schemes (de Jong et al. 1980; Leung et al. 
1984) or incorporating dust grains to the chemistry (Schilke et al. 2001). According to de 
Jong et al. (1980), the C^ ion, highly abundant in UV-processed cloud envelopes, would start 
the chain of reactions leading to CH2 and finally, to HCO. Confirming these expectations, 
HoUis & Churchwell (1983) found an empirical correlation between HCO and C"^ radio 
recombination line emission in two galactic clouds. Schilke et al. (2001) have proposed a 
model where HCO is the final by-product of the photodesorption or evaporation of solid 
formaldehyde on dust grain mantles (Westley et al. 1995). Although the values of X(HCO) 
predicted by Schilke et al. (2001) are still one order of magnitude below what is typically 
found in PDR, these results are encouraging. Based on a 1-dimensional model made for the 
Orion molecular cloud, Schilke et al. (2001) predict that the large photodissociation rate of 
HCO (van Dishoeck 1988) could be counterbalanced only for extinctions A^, >5-6. For a 
population of clouds uniformly bathed in the pervasive UV field of a starburst galaxy, the 
A^ limit should be risen by a factor of ~a few, making the equivalent condition on column 
density close to N(H2)>10^^cm~^. 

The HCO map of M82 provides direct evidence that the starburst event has heavily 
processed the bulk of molecular gas in this galaxy. The global fractional abundance of HCO 
(~410~^'^, averaged over ~650pc) can be accounted within a PDR scenario. Paradoxically, 
the spatial correlation between HCO-enriched clouds and HH regions is poor however. This 
is noticeable both at large scales (HCO/CO/HH nested ring morphology) but also at small 



- 7- 



scales (peaks of HCO/HII emission avoid each other). The poor correlation suggests that UV 
fields coming from the strongest HII regions of M82, lying in the inner ~400pc, have pho- 
todissociated the bulk of HCO in nearby molecular cloud envelopes. Previous observational 
evidence pointed out to a disrupted physical environment for the dense ISM in M82. Obser- 
vations of the atomic carbon [CI] ^^Pi and ^Pi ^'^Pq lines (Schilke et al. 1993; Stutzki 
et al. 1997) confirmed that the M82 [CI]/CO abundance ratio (~0.5) is higher than observed 
in non-starburst disks (it is ~.15 in our Galaxy). Moreover, the measured J=2-l/l-0 ratio 
can only be reconciled with clouds being small (with sizes ~lpc), only moderately dense 
(with densities ~10^-10^cm-=^), and hot (with temperatures ~50-100K) (Stutzki et al. 1997). 
Similar conclusions, based on Large Velocity Gradient (LVG) and PDR modeUing of ^^CO, 
^^CO, and C^^O fine ratios have been reported by Mao et al. (2000) and Weiss ct al. (2001). 
The bulk of molecular clouds in the inner ~400 pc of M 82 might have low densities (at some 
places in the disk as low as ^lO^cm"'^) and small sizes (~lpc). Assuming this scenario, 
the survival of HCO against photodissociation, requiring column densities N(H2)>10^^cm~^ 
might be difficult. Most remarkably, the location of HCO-enriched clouds in our map, iden- 
tified by their largest N{HCO)/N{H^^CO^) ratio (the outer ring and the central clump), 
coincides with Weiss et al.'s model predictions on where individual molecular clouds may 
have the largest column densities (where volume densities reach a few ~10^cm~^). Further- 
more, the intensity of UV fields is there comparatively lower (see Figure 3) . 

The reported differences between HCO and H^^CO+ maps in M82 cannot be easily 
explained with the present PDR models however, as both species are expected to be closely 
associated in PDR. Other unexplored scenarios cannot be excluded as plausible explanations 
for the enhancement of HCO: X-Ray or Cosmic Ray Dominated Regions chemistries might 
account for it in the case of M82 (Suchkov, Allen, & Heckman 1993). Moreover, we can 
not exclude the presence of dense molecular gas (with n(H2)>10^cm^^) in the central 1 kpc 
of M82 as already pointed out by Mao et al. (2000) and Weiss ct al. (2001), or shown by 
the detection of tracers of dense molecular gas (see HCN map of Brouillet & Schilke 1993). 
However, our results suggest that evaporation/destruction of molecular clouds by UV fields 
is highly efficient in the inner ~400pc of the galaxy. This efficiency shows spatial trends 
within the M82's disk but the reasons behind these variations remain to be understood. 

M 82 is the only galaxy where emission of HCO has been detected so far. HCO emission 
was also searched for in the starburst galaxy NGC253 by Garcia-Burillo et al. (2000). Their 
data, giving no detection, allowed to set a 3-a upper limit for the HCO/H^^CO+ ratio of 
< 0.12, i.e. a factor of 4-5 smaller than the value measured in the M82's ring. HCO is 
the only complex molecule showing a larger fractional abundance in M 82 than in NGC 253 
(Mauersberger & Henkel 1993). Furthermore, the remarkably different abundances and 
spatial distributions of SiO gas in NGC 253 and M 82 are suggestive of an evolutionary link 
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between these starbursts (Garcia-Burillo et al. 2000, 2001). The large HCO abundance 

in M82 may fit within this scenario, which considers the M82 starburst episode as more 
evolved than in NGC 253. The chemistry of molecular gas in NGC 253 is heavily influenced 
by the large-scale shocks and heating induced by a burst of pre-main sequence massive stars. 
This explains the high abundance of SiO and the significantly low abundance of HCO. The 
chemistry of molecular gas in M82 is dominated by the action of UV fields produced by 
more evolved massive stars giving rise to HII regions. The remnant of the starburst presents 
a PDR-hke chemistry favoring the presence of HCO, but forcing a low abundance for other 
complex molecules. 

This work received support from the Spanish DGES under grant number AYA, 2000- 
0927 and CICYT-PNIE under grant ESP2001-4519-PE. We acknowledge the IRAM staff for 
carrying the observations and help during data reduction. 

REFERENCES 

Achtermann, J. M., & Lacy, J. H. 1995, ApJ, 439, 163 

Alton, P. B., Davies, J. I., & Bianchi, S. 1999, A&A, 343, 51 

Bregman, J. N., Schulman, E., & Tomisaka, K. 1995, ApJ, 439, 155 

Brouillet, N., & Schilke, P. 1993, A&A, 277, 381 

de Jong, T., Boland, W. & Dalgarno, A. 1980, A&A, 91, 68 

Garcia-Burillo, S., Martin-Pintado, J., Fuente, A., & Neri, R. 2000, A&A, 355, 499 
Garcia-Burillo, S., Martin-Pintado, J., Fuente, A., & Neri, R. 2001, ApJ, 563, L27 
HoUis, J.M. & Churchwell, E. 1983, ApJ, 271, 170 
Joy, M., Lester, D. F., & Harvey, P. M. 1987, ApJ, 319, 314 

Larkin, J. E., Graham, J. R., Matthews, K., Soifer, B. T., Beckwith, S., Herbst, T. M., & 
Quillen, A. C. 1994, ApJ, 420, 159 

Leung, C. M., Herbst, E., & Huebner W. F. 1984, ApJS, 56, 231 

Mauersberger, R., & Henkel, C. 1993, Astronomische Gesselschaft, 6, 69 



-9- 



Mao, R. Q., Henkel, C, Schulz, A., Zielinsky, M., Mauersberger, R., Storzer, H., Wilson, T. 
L. & Gensheimer, R 2000, A&A, 358, 433 

Rieke, G. H., Lebofsky, M. J., Thompson, R. I., Low, F. J., & Tokunaga, A. T. 1980, ApJ, 
238, 24 

Sage, L. J., & Ziurys, L. M. 1995, ApJ, 447, 625 
Sakai, S., & Madore, B. F. 1999, ApJ, 526, 599 

Schenewerk, M. S., Snyder, L. E., & Hjalmarson, A. 1986, ApJ, 303, L71 

Schenewerk, M. S., Snyder, L. E., HoUis, J. M., Jewell, P. R., & Ziurys, L. M. 1988, ApJ, 
328, 785 

Schilke P., Carlstrom J. E., Keene, J., & Phillips, T. G. 1993, ApJ, 417, L67 

Schilke P., Pineau des Forets G., Walmsley C. M., & Martin-Pintado J. 2001, A&A, 372, 
291 

Seaquist, E. R., Carlstrom, J. E., Bryant, P. M., & Bell, M. B. 1996, ApJ, 465, 691 

Seaquist, E. R., & Clark, J. 2001, ApJ, 552, 133 

Shopbell, P. L., & Bland-Hawthorn, J. 1998, ApJ, 493, 129 

Snyder, L. E., HoUis, J. M., & Uhch, B. L. 1976, ApJ, 208, L91 

Snyder, L. E., Schenewerk, M. S., & HoUis, J. M. 1985, ApJ, 298, 360 

Stutzki, J., Graf, U. U., Haas, S., Honingh, C. E., Hottgenroth, D., Jacobs, K., Schieder, R., 
Simon, R., Staguhn, J., Winnewisser G., Martin, R. N., Peters, W. L., & McMuUin, 
J. P. 1997, ApJ, 477, L33 

Suchkov A., Allen R. J., & Heckman, T. M. 1993, ApJ, 413, 542 

van Dishoeck, E. F. 1988, Rate Coefficients in Astrochemistry, ed.T. J. Millar, & D. A. 
Williams, Kluwer, Astrophys. Space Science Library, 146, 49 

Weiss, A., Neininger, N., Hiittemeister, S., & Klein, U. 2001, A&A, 365, 571 

Westley, M. S., Baragiola, R. A., Johnson, R. E., & Baratta, G. A., 1995, Nature, 373, 405 

Wills, K., Redman, M. P., Muxlow, W. B., & Pedlar, A. 1999, MNRAS, 309, 395 



This preprint was prepared with the AAS macros v5.0. 



-10- 



Fig. 1. — The top panel (a) shows the velocity-integrated intensity map of HC0(F=2-1) in 

the central region of M82 (contour levels: -0.09, 0.144 to 0.350Jy beam~^km by steps 
of 0.024Jy beam^^km s^^; (7=0. 048 Jy beam~^km s^^ ). A filled ellipse on the lower right 
corner pictures the (5.9"x5.6") synthesized beam. The dashed line traces the galaxy major 
axis (PA=70°). (AX, AY) offsets, along the major and minor axis respectively, are referred 
to the dynamical center of the galaxy, marked by the central filled square. The bottom 
pannel (b) shows the spectra integrated over the squared regions covering the eastern and 
western lobes as shown. The hyperfine components of HCO detected are labeled by their 
"F-jump". H^3CO+(J=1-0) and SiO (v=0, J=2-l) spectra are also plotted. 

Fig. 2. — Top panel (a) shows the comparison of the integrated-intensity maps of HCO 
(F=2-l) (contours: 0.144 to 0.350Jy beam~^km s^^ by steps of 0.024Jy beam~^km s^^) 
and H^^CO+(J=1-0) (gray scale: 0.280-1 Jy beam~^km s~^). Bottom panel (b) shows 
comparison of HCO (F=2-l) and ^^C0(2-l) integrated-intensity maps (gray scale for CO: 
200-550Jy beam-^km s^^). 

Fig. 3.— Map of the HCO(F=2-1)/H^^CO+(J=1-0) integrated-intensity ratio (gray scale 
from 0.15 to 0.5) compared with the [Nell] flux map of Achtermann & Lacy (1995), in 
contours linearly scaled from 10% to 90% by steps of 10% of the peak value . 
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